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We dedicate this article to our teacher S.P. Novikov 
on the occasion of his 10th birthday 



In the seminal paper [21] in 1974 S.P. Novikov, in particular, established 
that the spectral curve of the one-dimensional periodic Schrodinger operator 

H = — — 7j + u(x) 
dx A 

is preserved when the real-valued potential u(x, t) evolves via the Korteweg- 
de Vries (KdV) equation and that for finite-zone (finite gap) potentials the 
classical conservation laws, i.e. the Kruskal-Miura integrals, are described in 
terms of branch points for this curve. The spectral curve T is a hyperelliptic 

A 2 = Q{E) 

where 

Q(E) = {E - Eq) . . . (E — E 2N ) 

is a polynomial of degree 2N + 1 for iV-zone potentials. It was proved in [21] 
that finite-zone potentials are exactly solutions of the Novikov equations, 
i.e., stationary points of higher KdV flows and their linear combinations, 
and that the KdV flow on the set of TV-zone potentials reduces to a com- 
pletely integrable finite-dimensional Hamiltonian system for which the ends 
of the stability zones, i.e., Eq, . . . , E 2 n, supply the necessary family of first 
integrals. 

The article [24] was the starting point for the development of the finite 
gap integration theory in which the spectral curves play the main role. 

In this article we consider the deformation of the spectral curve via the 
periodic equations of the Melnikov type and we show that although the 
spectral curve is not preserved it is deformed in such a manner that it still 
gives many conservation laws for the system. 

1. Introduction 

We recall that the KdV equation 

Ut = 6uu x — u xxx 
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has the Lax form 

(1) H t = [H,A] 

and as the spectral curve of H parameterizes the Bloch (— Floquet) functions 
which are formal eigenfunctions of H (here we do not mean that they lie in 
some nice functional space) and the monodromy operator Tf(x) = f(x + T) 
where T is the period of u(x): 

(2) H^j = EtP, fip{x) = e^ T ip(x) 

where [i is the quasimomentum which is defined on the spectral curve: fi = 
/i(A, E). The t-deformation of u results in the deformation of ip via the flow 

ip t = Atp. 

Another form of soliton equations instead of the Lax form is the Manakov 
triple: 

(3) H t = [H,A]+BH 

where A and B are differential operators. The main example is given by 
the Novikov-Veselov (NV) equations |27j for which H is a two-dimensional 
Schrodinger operator: H = dd + u. The spectral curve of H on the zero 
energy level T parameterizes only Floquet functions corresponding to the 
zero energy level: 

Hip = 0, 

(4) j>(x + T lt y) =e^ lTl V>(x,y), 

iP(x,y + T 2 ) =e i ^ T *iP(x,y) 

where T\ and T 2 are the periods of u. This curve was first introduced by 
Dubrovin, Krichever, and Novikov in [4j where the inverse problem at one 
energy level for two-dimensional Schrodinger operators was posed and solved 
for finite-zone operators (the spectral data for potential operators, i.e. with 
no magnetic field, were later distinguished in [26J). Therewith the Floquet 
functions are deformed again via ipt = Aip and hence the spectral curve is 
again preserved and may be considered itself as a conservation law. Another 
equation of such triple form is the modified Novikov-Veselov equation for 
which H is a two-dimensional Dirac operator and which being introduced 
by Bogdanov found applications in the surface theory |29j . 

Another generalizations of the Lax equations was proposed by Melnikov 
[18j and later was also derived by Kuznetsov and Zakharov [32J . The general 
form of these equations is the following extension of the Lax form: 

H t = [H,A] + C 

where 

N 
n=l 
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is the sum of differential operators Cj with coefficients depending on solu- 
tions (f>n, . . . , (fiiki of the auxiliary linear problems 

Hipi k = \iip ik , k = l,...,ki. 

Very frequently these equations are called the equations with self-consistent 
sources, each of them has a soliton predecessor of the form H t = [H, A] and, 
for example, the KdV equation with self-consistent force takes the form 

1 3 N 

(5) u t = -u xxx - -uu x + 2d x ^ipkix, t)ipl (x, t) 

k=i 

where ipk(x,t), ipl(x,t) are some solutions of the auxiliary linear problem 

{-dl + u)ip k = E k ip k , 

(-d 2 x + u)^ k = E k r k . 

To obtain a well-defined dynamics it is natural to assume that the products 
of the eigenfunctions in ([5]) are bounded. The simplest choice for the periodic 
problem is the following: ipk( x ) is a Bloch eigenfunction and ipt(x) is the 
Bloch eigenfunction with the inverse Bloch multipliers. 

The theory of such equations was developed in series of papers by Mel- 
nikov [191 1201 12D 1211 [23] and others mostly for the case of functions fast 
decaying at infinity. 

In this article we show that in difference with soliton equations the spec- 
tral curve is not preserved by these systems however it still gives many 
conservation laws. 

2. The spectral curve 

The systems ([2]) and ([3]) do not have solutions for all possible values of 
constants, i.e. for all E and \x in the former case and for all \i\,\X2 in the 
latter case. In fact, such solutions exist if and only if these constants satisfy 
some analytical condition ("the dispersion laws"): 

F(E,v) = 0, G(/x 1 ,/x 2 ) = 0. 

Each equation describes a complex curve T C C 2 and to each point of T 
there corresponds a linear space of solutions to the corresponding equation, 
([21) or @ . This picture was drawn in physical terms in [25] and two different 
methods for the justification of it were proposed by Krichever and the second 
author (I.A.T.) (see [13 ED [8] ) . 

Now to obtain the spectral curve we have to consider the ^-bundles 
formed by solutions to ([2]) or @ and normalize T at such a manner that the 
pull-back of the ^-bundle onto under the projection 

form a bundle with fibers of constant dimension. We refer for details to |31] 
and here demonstrate this procedure by an important original example. 
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Example. [24J For the one-dimensional Schrodinger operator with a 
real-valued potential the multipliers of T are defined on a Riemann surface 
(a complex curve) T 

A 2 = Q{E) 

where Q(E) is an entire function with infinitely many zeroes. All zeroes 
lie on the real line. To every point P = (E, A) G T where Q / there 
corresponds a one-dimensional space of solutions to ([2]). Let fi'gi satisfy 
the following conditions 

(1) Q has a zero at E' of multiplicity two; 

(2) to the point (E',0) there corresponds a two-dimensional space of 
solutions to ([2]). (This, in particular, implies that this is a double 
point on r.) 

Let us unglue this double point and obtain another Riemann surface T'. 
Then the '(/'-bundle over T is pulled back to a bundle ip 1 over T' with one- 
dimensional fibers at the preimages of {E 1 , 0). Moreover this bundle is holo- 
morphic near these points. We have 

• if all zeroes of Q except finitely many satisfy conditions 1 and 2 
above, then after ungluing all corresponding double points we obtain 
a Riemann surface of finite genus and the one-dimensional ip- 
bundle over it. The surface is defined by the equation 

A 2 = Q{E) 

where Q is a polynomial of odd degree, say 2N + 1. 

It is said that this operator is finite-zone (or finite gap), and it all 
zeroes of Q are simple it is said that it has N zones (gaps). There 
is a function ip{P, x) meromorphic in P € with the following 
asymptotic 

(6) V ~ ^ X as E -» °°- 

Therewith the complex curve is compactified to an algebraic 
curve by adding the point E = oo and tp becomes a meromorphic 
function on T^p \ {E = 00} with the essential singularity (|6|) at 
E = 00. 

Here we remark that itself may have singularities and, in fact, there 
is a tower of projections 

Lnorm ^ T^, > T 

where r norm is the normalization of T. The multiplier mapping which cor- 
responds to a point the set of "multipliers" : 

M:T^C 2 , M(E,X) = (E,fi) 

is naturally ascends to this tower. 

Above we explain how the spectral curve arises from the spectral theory 
of differential operators. 
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However the strongest method for constructing exact periodic (and also 
quasi-periodic) solutions of solitons equation, i.e. the Baker- Akhiezer func- 
tion method [13], starts with an introduction of an algebraic curve T and of 
a function ip (which may be a vector or even matrix function) with asymp- 
totics of the kind of ([6]) at several points of T. It is assumed that tp is 
defined by some additional data uniquely. The function ip is a formal eigen- 
function of some operator H which is uniquely reconstructed in terms of of 
algebraic functions corresponding to V from tp. Such a function ip is called 
the Baker- Akhiezer function of H , and the soliton dynamics ([I]) or ([3]) is 
linearized in terms some data coming in the definition of tp and this leads 
to explicit algebra-geometrical formulas for so-called finite gap solutions of 
soliton equations. The spectral curve T is preserved by the flow, i.e. the 
flow is isospectral. 

Therewith for operators H with periodic coefficients and with a nice spec- 
tral theory (i.e. for which the existence of the dispersion laws may be es- 
tablished) r = r^, and tp is the section of the -0-bundle. 

In this article we show that 

• in contrast with soliton equations the periodic equations of the Mel- 
nikov type may be almost isospectral, i.e. it may preserve A4(T^) 
and deform T^. 

The first example of this effect was found by us in [8] and we expose it in 
the next section. 

3. The conformal flow for the (Weierstrass) potentials of tori in 

M 3 and M 4 

The author's interest to the study of the Melnikov-type equations was 
partially motivated by the problem of conformal invariance of the higher 
Willmore functionals. 

By the generalized Weierstrass method, any torus in IR 3 is described in 
terms of the zero-eigenfunction tp: 

Vtp = 0, 

of a two-dimensional periodic operator 

V=( °- 9 \ + ( U °) 

where the potential U is real-valued and any torus in M 3 is described in 
terms of two solutions (p, ip to the equations 

Vtp = 0, V v cp = 

where 

\ -d J + V U J ' V -d J + V U J 

are to conjugate periodic operators (see, for instance, [31]). The spectral 
curve of V is naturally defined (see (JH and §2) and contains in itself 
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the information of the Willmore functional which is defined for all closed 
surfaces immersed in M 4 as follows 

W(M) = [ |H| 2 d/i 
Jm 

where H is the mean curvature vector and d\x is the induced volume. 

This functional is invariant with respect to conformal transformations of 
the ambient space, i.e. if we have a conformal transformation / : IR 4 — > R 4 
which maps a compact surface without boundary M into a compact surface, 
then 

W(M) =W(/(M)). 

This follows from the conformal invariance of the form (|H| 2 — K)dfi where K 
is the Gaussian curvature and the Gauss-Bonnet theorem by which j Kd[i 
equals 2irx(M), i.e. the topological quantity. 

The soliton local deformations of surfaces in M 3 and M. via the modified 
Novikov-Veselov (mNV) equation and the Davey-Stewartson (DS) equation 
were introduced by Konopelchenko OHO]- It appears that they preserve the 
tori globally and therewith preserve the Willmore functional as well as the 
spectral curve \29\ 131]. hence it is natural to treat higher conservation laws 
of these hierarchies as higher Willmore functionals. 

The conformal invariance of the Willmore functional led the second author 
(I.A.T.) to the conjecture that these higher Willmore functionals and the 
spectral curve for tori in M. 3 themselves are conformally invariant [30J. 

It was rather soon established by the first author (P.G.G.) and M.U. 
Schmidt [7] who considered the conformal flow, i.e. the Melnikov type flow, 
induced on the potential U by continuous conformal transformations: 

u T = KM 2 -KM 2 

where the torus is defined via the Weierstrass formulas by ifi = (V , i,V ; 2) T - 
Under this deformation the "^-function on the spectral curve evolves in such 
a manner that the quasimomenta are preserved. 

In [8] we analyzed carefully this situation for the more general case of tori 
in R 4 . It appears that the conformal flow on U which corresponds to the 
following generator of the conformal group 

df x — .2./ x . d-j-x — x . 

d T x 3 = (x 3 ) 2 - (x 1 ) 2 - (x 2 ) 2 - (x 4 ) 2 , 

d r x 4 = 2x 4 x 3 

has the Melnikov form: 

d T U = (fii)i - if2lp2, 

(7) - 

8 T U = (flip! - (£2^2 

where ifi = (ipi,ip2) T an d tp = {(fi,^)^ define a torus in M 4 via the gener- 
alized Weierstrass formulas. It appears to be isospectral in the sense that 
all multipliers are preserved. However we knew about several explicitly 
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computed examples of the Weierstrass representations of tori which are the 
Clifford torus in S 3 : 



and its stereographic projection into M 3 [31J. In these cases the spectral 
curves I" 1 ,/, are different: the complex projective line CP 1 in the former case 
and CP 1 with two pairs of points glued into two double points. However 
both tori are connected by a continuous conformal transformation of M 4 . A 
detailed analysis led us to the following conclusion: 

• the conformal flow (|7|), i.e. a particular case of Melnikov deforma- 
tions of periodic operators, is only almost isospectral, i.e. preserve 
the multipliers — the complex curve A4(T^) — and deform the spec- 
tral curve r^, . In this particular case the deformation of I\a consists 
in gluing and ungluing double points. 

• since the higher integrals of the mNV and the DS hierarchies are 
described in terms of M(T^j), these integrals are preserved and give 
us spectral conservation laws of the conformal flow. 

4. The Baker Akhiezer function and kernel and the (ip,ip*)- 
representation of equations 

Let us recall the definition of the Baker-Akhiezer function for the KP 
equation [15] . 

Let r be a smooth Riemann surface of genus g with the following data: 

(1) a divisor of poles D = 7x + . . . + 7 g ; 

(2) a distinguished point P with a local parameter z = 1/A. 

The Baker-Akhiezer function "0(7) t) depends on the spectral parameter 
7 € T and of infinite set of real variables x = t±, y = t2, t = ts, £4, £5, . . . , 
t = (x, y, t, £4, £5, . . .). To avoid analytic problems it is convenient to assume 
that t has only finite number of nonzero entries. 

For generic t there exists an unique function of 7 G T such, that: 

(1) V(7> t) i s meromorphic in 7 outside P with simple poles at 71,. . . ,7 g . 



If u is periodic in x and y, then T is the spectral curve (on the zero energy 
level) of the operator d y — d% + u(x, y) [15] . 






(9) -i> XX (\,t)+ijy(\,i)+u(t)ij(\,t) =0. 
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Let us assume that T is a hyperelliptic surface such that A 2 is a global 
meromorphic function on T with exactly one second-order pole at P. Then 



exp 



k>0 



4>(i,x,t,t 5 ,t 7 , . . .) 



and t/j(j,t), t = (x,t,t§,t7, . . .) is the Baker-Akhiezer function of the KdV 
hierarchy [5j. We shall omit the tilde sign in the KdV formulas. In the KdV 
case we have 

(10) - ip xx (X, *) + «(t)^(A, i) = -A^(A, t) 
instead of (|9]). 

In [3] Cherednik has shown that all flows from the KdV hierarchy are 
obtained as the expansion coefficients in A^ 1 near A -1 = for the following 
A-dependent nonlocal equation: 

(11) u T = 2d x (if) k (\x)ip k (-\,x)). 
Here we assume that all times except x are equal to 0. 

Theorem 1. Let the source functions ip k and ipt in the right-hand side of 
([5]) be the restrictions of the Baker-Akhiezer function at some points ofT: 

Then ([5]) can be represented as the following linear combination of the flows 

C2P- 



2SU 



-res I (A 3 ^(A, x, r)ip(— A, x, r)d\)+ 

N 

+ ^2 1p{^k,X, T)lp(-\ k ,X, T 
k=l 

All the higher KdV flows are isospectral and form a commutative algebra. 
Typically the complete algebra of symmetries for soliton equations is non- 
commutative and contains both isospectral and non-isospectral flows (see 
|28| for further references). 

Orlov and Schulman suggested a generic approach for studying the sym- 
metry algebra based on the so-called infinitesimal dressing [28]. It particular, 
in [28] it was shown that generators K mn [u] of the algebra of all KdV and 
KP symmetries is obtained by expanding the flow 

(12) u T = 2d*(V>fc(A,i)^GMl) 

near the diagonal A = fi: 

( i \ m / 1 r ~ 



A 



/i 



at the point P where A = fx = oo . For n = we have the standard KP 
(KdV) hierarchy. The n = 1 coefficients generate the Virasoro algebra of 
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non-sospectral symmetries, the n > 2 symmetries are not compatible with 
the KdV reduction. 

Here t) is the wave function for all auxiliary linear operators of the 
KP (KdV) hierarchy, t = (x = t\,t = £3,^5, . . .) or t = (x = t\,y = t2,t = 
£3, £4, . . .) denotes the full set of KdV (KP) times, and ip*(\, t) satisfy the 
formal conjugate linear problems. In the KdV case all auxiliary problems 
are self-adjoint, therefore 

V*(A,t)=V(-A,t). 

An arbitrary source function may be expanded in terms of eigenfunctions 
products. Such expansions play a critical role in the perturbations theory for 
soliton equations. The periodic perturbation theory for 1-dimensional finite- 
gap potentials and for the 2-dimensional finite-gap at one energy potentials 
was developed by Krichever [HI [15]. In particular, he pointed out that 
it is natural to treat the conjugate Baker-Akhiezer function ip*(\,t) as a 
holomorphic 1-form in the spectral parameter 7 on T\P. It is defined by 
the following analytic properties: 

(1) ifj*('j,i) is an 1-form in 7, i.e. in local coordinates it reads as 
ip* (A, t) = </>* (A, i)d\, where ^* (A, t) is an analytic function. 

(2) tp* (7, t) is holomorphic in 7 outside P with simple zeroes at 71 ,. . . ,j g . 



(3) ril,t) =exp 



(1 + o{l))d\ as 7 ~ P. 



fc>0 

The action of the Virasoro algebra symmetries on the finite-gap KP solu- 
tions (these symmetries generically result in non-trivial deformations of the 
complex structures on the spectral curves) was studied by Orlov and the first 
author in [6]. In particular, in [6] it was shown, that the infinitesimal de- 
formations of the Baker-Akhiezer function corresponding to the generators 
(|12p (infinitesimal Darboux transformations of the finite-gap KP solutions) 
are naturally written in terms of the so-called Cauchy-Baker-Akhiezer kernel 
uj(X, /J,,t): 

The kernel u>(X,fj,,t) is defined by the following analytic properties: 

(1) uj(\, is a meromorphic function in A and a meromorphic 1-form 
in /i on T\P. 

(2) For a fixed \i the function u>(\, fi, t) has simple poles at the points fx 
71, •••,7s- 

(3) For a fixed A the 1-form u>(\, fx, t) has simple zeroes at 71,. . . ,j g and 
a simple pole at A. 

(4) Lu(X,fi,i) = -^j + 0(1) near the diagonal A = fi. 

(5) For a fixed /j, the function cj(A, ^Aexp 
at the point A = P. 



k>0 



is regular in A 
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(6) For a fixed A the 1-form u(X, fx, rj/i 1 exp 



,fc>o 



is regular in fi 



at the point ji = P. 

For t = this kernel coincides with the Cauchy kernel on Riemann surfaces 
used by Koppelman [Hj. For data generating regular potentials u(i) the 
following explicit formula was suggested in [B] : 

u(\,fj,,x,y,t 3 ,t4, ■■■) = 



V>(A, x', y, t 3 ,U, . . .)ip*{fi, x', y, t 3 , t 4 , . . .)dx' . 



The upper limit of the integral depends on the quasimomenta at the points 
A, jj, and is chosen to make the integral convergent. An analogous represen- 
tation for the Cauchy kernels on Riemann surfaces for systems with discrete 
x was suggested earlier by I.M.Krichever and S.P.Novikov in [16j. 

5. The periodic Kadomstev Petviashvili equation with a self- 
consistent source 

To integrate the KP equation with the self-consistent sources we have to 
consider spectral curves with additional double points. Such curves corre- 
spond to the solitons on the finite-gap background [12]. We assume that 
we have the same spectral data as in the Section 4 plus 2N marked points 
Bh, R_, k = 1, . . . ,N. Denote the local parameters near these points by 
A. The Baker- Akhiezer depends on N extra real parameters ri, . . . , r^v, 
t — . . . ,tjv) and has the following analytic properties: 

(1) ipipf, t, r) is meromorphic in 7 outside P with g + N simple poles at 
71,. . . ,7 9 , R\,. . . , R%. 

(2) res (A, t, f)d\ = r k ^(R k _,t, f). 

\=R I 1 



(3) ^(7,f,r) 



exp 



'l + o(l)) as 7 ~ P. 



Lfc>0 

The properties of the conjugate Baker-Akhiezer 1-form are the following: 
(1) ip*(j,t,f) is meromorphic in 7 outside P with simple zeroes at 
71,. . . ,7 S and simple poles at RL,. . . , R^. 



res 



**(A,t,f) 

A=_R fe 



(3) rilXr) 



exp 



\=R k + 

(l + o(l))dA as7~P. 

The corresponding potential u{t, t) is defined by the formula ([8]) 

braic genus g with the 



fe>0 



Theorem 2. Let F be a Riemann surface of 
following KP data: 

(1) a divisor of poles D = 71 + . . . + 7^; 

(2) a distinguished point P with a local parameter z = 1/A; 
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(3) an additional collection of2N points R\, R^_, k = 1, . . . , N. Denote 
the local parameters near R_ by X. 

Then potential u(t, f) defined above satisfy the following equations with self- 
consistent sources: 



(13) 



8u(t, t ) _ 2Q <KRt , t, r)r (A, t, t) 



dX 



Proof. The Cauchy-Baker-Akhiezer kernel uj(X, fi,t,r) corresponding 
to this spectral data has the following analytic properties: 

(1) fi,t,r) is a meromorphic function in A and a meromorphic 1- 
form in /i on T\P. 

(2) For a fixed \i the function cj(A, ^, t, f) has simple poles at the points 
li, 7i,. . . ,7 9 , R\,.^ . ,R*j[ . 

(3) res uj(X, fj,,t,r)d\ = Tu>(R k _, fi,t,r). 

\=R* 

(4) For a fixed A the 1-form w(A,/i,t,r) has simple zeroes at 71,. . . ,7 S 
and simple poles at A, . . , R_ . 

(5) res lj(X, fi,t,r) = —TkLd(X,t,T)/dii 
n=R k _ 



H=R% 

(6) lu(X, ji, t, t) = -jj^j + 0(1) near the diagonal A = /i. 



E **** 

k>0 



E » k th 

k>0 



is regular in 
is regular in 



(7) For a fixed \x the function u;(A, fi, t, r)A exp 
A at the point A = P. 

(8) For a fixed A the 1-form cj(A, fi, t, t)/x _1 exp 

/x at the point /i = P. 
It is easy to check (see[B]): 

(14) d x uj(X,fi,t,T) = -ip(X,t,r)ip*(n,t,T), 

d y uj(X, nX?) = ~i>x{\ t, r)xp*(fi, t, t) + ip(X, t, r)^*(/x, t, r). 
We use the following formula: 

u(X,n,t,T) 



(15) 



d Tk tp(X,t,T) 



dfi 



^(Rl,t,r). 



To prove (|15p is sufficient to check that the right-hand side has the correct 
analytic properties. In particular, the relation 

^(A, t, r)dX = T k d Th i){R k _ , t, t) + i){R k _ , t, t) 

X=R I 1 



d Tk res 



follows from the following expansion near the point A ~ R+: 



U>(\,fJ,,t,T) 

dfi 



_uj{R k _,nXT) 
dfi 



ti=R* 



X-Ri) 



+ 0(1). 
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This expansion can be easily derived from the properties [3] and [6] of the 
Cauchy-Baker-Akhiezer kernel. 
Taking into account, that 

(16) ^xx(A)-^ y (A)-«V(A)=0, VL(M)+^(M)-^*(M) = 

we obtain the following formula for the variation of u = u(t, f) 



(17) 



Vw fe (A) - tp yTk (A) - mp Tk (A) 



Substituting ([15]). CO} into dTTJ) we obtain ([13]). We see, that the right-hand 
side of (|17p turns out to be A-independent, therefore the deformation (|15p 
of the Bloch function is admissible. 

Corollary 1. Denote by u(x,y,t), ^(7, x, y, t), i>* (7 , x , y , r) £/ie functions, 
obtained from u(t, f), ij}{^,t,r), ip*(^,t,T) by the following linear substitu- 
tion: 

(18) ti = x, t 2 = y, t k = c k r, k = 3, . . . , M, t fc = 0, k > M, 

Tk = a k + far, k = l,...,N. 

Then u(x, y, t) solves the following Melnikov-type equation: 

du{x,y,r) 



(19) M N 



dr 

$(Rt,x,y,T)i)*(\,x,y,T) 



dX 

k=3 k=l 



where K k [u] denotes the k-th flow from the standard KP hierarchy, the func- 
tions in the right-hand side satisfy klb]) 

Remark. It is easy to notice from the previous formulas that0 
if at t = tiie function u is periodic in x and y, tp is the Floquet eigen- 
function of the operator L = d y — d\ + u, and for k = 1, . . . , N the products 
ip(R_, x, y, t)^*(R^_, x, y, r) are periodic, then the evolution [W\) preserves 
the periodicity of u and the multipliers of tp. 
We conclude that 

• In fact, the Baker- Akhiezer function in Theorem [2] is defined on the 
spectral curve with double points. The double point (R+,Rt) is 
unglued if and only if r k = 0. For generic r all double are present, 
i.e. , the spectral curve is obtained from T by pair-wise gluing points 
(R\,Rl),..., (R%, RN), and the fc-th double point is unglued when 
a k + flk T = 0. If the initial spectral curve is regular for r = 0, then 
equations with self-consistent sources immediately generate double 
points, which remains unglued for almost all times. 



^The derivation of the similar fact for the conformal flow (see Section 3) is exposed in 
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Another example of such an effect is given by the conformal flow 
(see Section 3). 

This observation demonstrates the principal difference between the Mel- 
nikov-type equations and the standard hierarchies like KdV, nonliear Schro- 
dinger (NLS), sine-Gordon, or KP. The standard hierarchies are isospectral, 
therefore the evolution can not generate double points or unglue the existing 
double points in a finite time. 

In the case of the self-focusing NLS equation the spectral curves with 
double points may correspond to regular space-periodic solutions, associated 
with so-called whiskered tori. All double points for such solutions remain 
glued for all values of t but become unglued in the limit t — > ±00. In [TJ [2] 
it was shown that for a periodic solution corresponding to a smooth curve 
the generation of double points after arbitrary small perturbations results 
in numerical chaos. 

6. An annihilation of a soliton in a finite time 

Let us discuss the simplest example - the one-soliton solution of the KdV 
equation. 

The wave function tp(X, x, c) of a one-dimensional Schrodinger operator: 
—ip" + uip = Eip, u = u(x,c), E = (iX) 2 , 
has the following form 

^(X,x,c)=e Xx (l + X ^ C ^ 



X + K 

The spectral curve T is the Riemann sphere with a double point: A = — k 
and X = k are glued together. We assume that the divisor point is located 
at A = — k, therefore we have: 

les ijj(X,x,c) = — cip(K, x, c), 

(this relation coincides with the property ([2]) of the Baker-Akhiezer function 
from the Section 5) and 

, , -2ck e KX 
X{x,c) 



2ne~ KX + ce KX 
By (JHJ) we obtain 

— 16ck 3 

(20) u(x, c) = 2d x x{x, c) = -2 . 

[ZKe KX + ce KX ) 

The c-dynamics corresponds to the following choice of Baker-Akhiezer func- 
tion solutions: 

tl>i(x, c) = ib(n, x, c) = . 

' > > ! ) 2ne- KX + ce KX 

The conjugate Baker-Akhiezer function is defined by: 

Tp* (A, x, c) = A, x, c)dk. 
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A simple straightforward calculation shows, that 
(21) d c u(x,c) = —2d x ip 2 (K,x,c). 



Formula (I20p generates regular solitons for c > 0, singular solitons for c < 
and zero solution for c = 0. If c — > 0, the position of the soliton goes to the 
+oo. 

The standard KdV evolution of a soliton is given by (|2ip where c = c(t) 
is governed by: 



We see that 

• starting with sufficiently small c we reach the point c = at a finite 
time. At this moment the double point on the spectral curve vanishes 
and the soliton annihilates. 

By inverting the direction of evolution we obtain examples of of such 
effects as 

• a creation of a soliton (there is no soliton at c = and it exists as 
soon as c > 0); 

• a capture of a soliton (c — * k~ 3 as t — > oo), 

first observed by Melnikov |19j . 
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d tC (t) = K 3 c(t). 

Let us consider the following Melnikov-type flow: 




The corresponding evolution of c(t) is given by 

dtc = k 3 c — 1. 
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